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Abstract. Primordial nucleosynthesis, or Big Bang Nucleosynthesis (BBN), is one of the 
three evidences for the Big-Bang model, together with the expansion of the Universe and 
the Cosmic Microwave Background. There is a good global agreement over a range of nine 
orders of magnitude between abundances of 4 He, D, 3 He and 7 Li deduced from observations, 
and calculated in primordial nucleosynthesis. This comparison was used to determine the 
baryonic density of the Universe. For this purpose, it is now superseded by the analysis of 
the Cosmic Microwave Background (CMB) radiation anisotropics. However, there remain, a yet 
unexplained, discrepancy of a factor 3-5, between the calculated and observed lithium primordial 
abundances, that has not been reduced, neither by recent nuclear physics experiments, nor by 
new observations. We review here the nuclear physics aspects of BBN for the production of 
4 He, D, 3 He and 7 Li, but also 6 Li, 9 Be, n B and up to CNO isotopes. These are, for instance, 
important for the initial composition of the matter at the origin of the first stars. Big-Bang 
nucleosynthesis, that has been used, to first constrain the baryonic density, and the number 
of neutrino families, remains, a valuable tool to probe the physics of the early Universe, like 
variation of "constants" or alternative theories of gravity. 



1. Introduction 

There are presently three evidences for the Big-Bang Model : the universal expansion, the 
Cosmic Microwave Background (CMB) radiation and Primordial or Big-Bang Nucleosynthesis 
(BBN). The third evidence for a hot Big-Bang comes, indeed, from the primordial abundances 
of the "light elements": 4 He, D, 3 He and 7 Li. They were produced during the first ~20 minutes 
of the Universe when it was dense and hot enough for nuclear reactions to take place. These 
primordial abundances are compared to astronomical observations in primitive astrophysical 
sites. It is worth reminding that Big Bang Nucleosynthesis has been essential in the past, to 
first estimate the baryonic density of the Universe, p& = (1 — 3) x 10~ 31 g/cm 3 [T], and give an 
upper limit on the number neutrino families N v < 3 [2], both in the seventies. The number of 
light neutrino families is now known from the measurement of the Z° width by LEP experiments 
at CERN: N v = 2.9840±0.0082 [3j. The nuclear reaction rates have all been measured in nuclear 
physics laboratories or can be calculated from the standard theory of weak interactions. In that 
case, they are normalized to the experimental value for the lifetime of the neutron. Its precise 
value is still a matter of debate [3] T n = 880-884 s, but its uncertainty has only marginal effect on 
BBN. The last parameter to have been independently determined is the precise value of baryonic 
density of the Universe, which is now deduced from the observations of the anisotropies of the 
CMB radiation. It is usual to introduce rj, the number of photons per baryon which remains 



constant during the expansion, and is directly related to fib by f2b - ^ 2 =3.65x 10 r\ with 



n h -h 2 = 0.02249 ± 0.00062 and fi b = 0.04455 ± 0.0027 



(1) 



("WMAP only Seven Year Mean" [5]). The parameter h represents the Hubble constant (Hq) 
in units of 100 km/s/Mpc (i.e. h= 0.704 |5|) and fib= Pb/po,c the baryonic density relative to 
the critical density, po,c, which corresponds to a flat (i.e. Euclidean) space. It is given by : 



where G is the gravitational constant. It corresponds to a density of a few hydrogen atoms per 
cubic meter or one typical galaxy per cubic megaparsec (Mpc). This results (Eqs. HH2]) in a 
baryonic density which is just slightly above the range provided by Wagoner [1] in 1973! 

Hence, the number of free parameters in Standard Big Bang Nucleosynthesis has now been 
reduced to zero, and the calculated primordial abundances are in principle only affected by 
the moderate uncertainties in some nuclear cross-sections. It may appears that Big Bang 
Nucleosynthesis studies are now useless, but this is certainly not the case. First, even though 
the agreement with observations is good or very good for 4 He, 3 He and D, there is a tantalizing 
discrepancy for 7 Li that has not yet found a consensual explanation. Second, when we look back 
in time, it is the ultimate process for which, a priori, we know all the physics involved. Hence, 
departure from its predictions could provide hints or constraints on new physics or astrophysics. 

Besides the 4 He, D, 3 He and 7 Li isotopes, some minute traces of 6 Li, 9 Be, n B and CNO 
are produced by BBN. Observations of 6 Li in a few halo stars have renewed the interest for 
this isotope and the nuclear uncertainties concerning its production. The evolution of the first 
generation (Population III) of stars could be influenced by the amount of primordial CNO 
elements, as hydrogen burning can proceed either through the slow p-p chain, or through the 
more efficient CNO cycle. We will hence review the nuclear aspects of the primordial production 
of element up to oxygen. 

2. The cosmological elements 4 He, D, 3 He and 7 Li 

2.1. Abundances of the cosmological elements 

During the evolution of the Galaxy, complex nucleosynthesis takes place, mainly in massive 
stars which release matter enriched in heavy elements into the interstellar medium when they 
explode as supernovae. Accordingly, the abundance of heavy elements in the gas, at the origin 
of star formation, increases with time. The observed abundance of metals (in astrophysics, the 
chemical elements beyond helium) is hence an indication of its age: the oldest stars have the 
lowest metallicity. To derive the most primitive abundances one has first, to extract them from 
observations of astrophysical sites which are thought to be non evolved and second, extrapolate 
them to zero metallicity. 

Primordial lithium abundance is deduced from observations of low metallicity stars in the 
halo of our Galaxy where the lithium abundance is almost independent of metallicity, displaying 
a plateau [BJ. This constant Li abundance is interpreted as corresponding to the BBN 7 Li 
yield. This interpretation assumes that lithium has not been depleted at the surface of these 
stars so that the presently observed abundance is supposed to be equal to the initial one. The 
small scatter of values around the "Spite plateau" is an indication that depletion may not have 
been very effective. Astronomical observations of metal poor halo stars have led to a relative 
primordial abundance [7] of: 



Po,c = 



8ttG 



1.88 h 2 x 10~ 29 g/cm 3 or 2.9 h 2 x 10 11 M /Mpc 3 



(2) 




(3) 



Note also that observationally challenging detections of 6 Li has been reported [8] to a level of 
~ 1CP 2 below the Spite plateau value. The presence of a 6 Li plateau is however not confirmed as 
only few (2-3) stars seem to present significant 6 Li abundances on their surfaces [9| . For a recent 
review of the latest Li observations and their different astrophysical aspects, see Refs. |10|. Illj. 
and the proceedings of the 2012 Workshop Lithium in the Cosmo^. 

Contrary to 7 Li which can be both produced (spallation, asymptotic giant branch (AGB) 
stars, novae) and destroyed (in the interior of stars), deuterium is a very fragile isotope, that can 
only be destroyed after BBN. Its most primitive abundance is determined from the observation of 
absorption lines in clouds at high redshift, on the line of sight of more distant quasars. Very few 
observations of these cosmological clouds are available and a weighted mean |12j (and references 
therein) of this data yields a D/H abundance of: 

D/H = (3.02 ±0.23) x 1(T 5 . (4) 

After BBN, 4 He is produced by stars. Its primitive abundance is deduced from observations 
in HII (ionized hydrogen) regions of compact blue galaxies. Galaxies are thought to be formed 
by the agglomeration of such dwarf galaxies which are hence considered as more primitive. 
The primordial 4 He abundance Y p ( 4 He mass fraction) is given by the extrapolation to zero 
metallicity but is affected by systematic uncertainties such as plasma temperature or stellar 
absorption. Using the data compiled in Ref. [13], it was found [14^ 1 15 j that: 

Y p = 0.2534 ± 0.0083. (5) 

Contrary to 4 He, 3 He is both produced and destroyed in stars so that the evolution of its 
abundance as a function of time is not well known, and has only been observed in our Galaxy 

3 He/H = (1.1 ± 0.2) x 10~ 5 . (6) 
Consequently, comparison with 3 He abundance from BBN is subject to caution |17| . 

2.2. Nuclear reactions for 4 He, D, 3 He and 7 Li nucleosynthesis 

Unlike in other sectors of nuclear astrophysics, nuclear cross sections have usually been directly 
measured at BBN energies (~100 keV). There are 12 nuclear reactions responsible for the 
production of 4 He, D, 3 He and 7 Li in Standard BBN. There are many other reactions connecting 
these isotopes, but their cross sections are too small and/or reactants too scarce to have any 
significant effect. Even among these 12 reactions, a few of them are now irrelevant (see below) 
at WMAP baryonic density. 

The weak reactions involved in nf>p equilibrium are an exception; their rates [IB] come from 
the standard theory of the weak interaction, normalized to the experimental neutron lifetime. 
Until very recently, the averaged value of 885.7±0.8 s was recommended by the Particle Data 
Group, but new measurement lead to significantly lower values. While it has not yet been 
possible to solve this discrepancy [I], a reevaluation of the recommended value: 880.1±1.1 s has 
been proposed [TU], awaiting experimental clarification. 

The 1 H(n,7) 2 H cross section is also obtained from theory |20j but in the framework of Effective 
Field Theory. 

For the ten remaining reactions, 2 H(p,7) 3 He, 2 H(d,n) 3 He, 2 H(d,p) 3 H, 3 H(d,n) 4 He, 
3 H(a, 7 ) 7 Li, 3 He(d,p) 4 He, 3 He(n,p) 3 H, 3 He(a, 7 ) 7 Be, 7 Li(p,a) 4 He and 7 Be(n,p) 7 Li, the cross 
sections have been measured in the laboratory at the relevant energies. Formerly, we used 



1 http://www.iap.fr/lithiuminthecosmos2012/index.html, proceedings will appear in Memorie della Societa 
Astronomica Italiana Supplementi. 



Table 1. Abundance sensitivity to reaction rates: d logy/5 log < av > at WMAP baryonic 
density. 



Reaction 


4 He 


D 


3 He 


7 Li 


E Q { AE /2) 












(MeV) 


n^P (r" 1 ) 


-0.73 


0.42 


0.15 


0.40 




^(n^H 





-0.20 


0.08 


1.33 




2 H(p, 7 ) 3 He 





-0.32 


0.37 


0.57 


0.11(0.11) 


2 H(d,n) 3 He 





-0.54 


0.21 


0.69 


0.12(0.12) 


2 H(d,p) 3 H 





-0.46 


-0.26 


0.05 


0.12(0.12) 


3 H(d,n) 4 He 








-0.01 


-0.02 


0.13(0.12) 


3 H(a, 7 ) 7 Li 











0.03 


0.23(0.17) 


3 He(n,p) 3 H 





0.02 


-0.17 


-0.27 




3 He(d,p) 4 He 





0.01 


-0.75 


-0.75 


0.21(0.15) 


3 He(a,7) 7 Be 











0.97 


0.37(0.21) 


7 Li(p,a) 4 He 











-0.05 


0.24(0.17) 


7 Be(n,p) 7 Li 











-0.71 





the reaction rates from the the evaluation performed by Descouvemont et al. [2T]. However, 
more recent experiments and analysis have lead to improved reaction rates for several important 
reactions. 

To point out the most important reactions, we display in Table [1] the sensitivity of the 
calculated abundances (Yi with i = 4 He, D, 3 He and 7 Li) w.r.t. to a change in the 12 
reaction rates by a constant factor. We define the sensitivity as d logy/5 log < av >, based on 
the assumption that the nuclear cross section uncertainties are now dominated by systematic 
uncertainties that affect their normalization rather than by statistics. These values were 
obtained, at the WMAP baryonic density, by a parabolic fit of Ay /y for ±15% variations of the 
reaction rates. The last column represents the Gamow window at BBN typical temperatures. 

This table can be used as a guide for further experimental efforts. We see for instance that 
at WMAP baryonic density, the 3 H(a, 7) 7 Li and 7 Li(p,a) 4 He reactions play a negligible role. 
At WMAP baryonic density, 7 Li is produced indirectly by 3 He(a, 7) 7 Be, that will, much later 
decay to 7 Li. 

The sensitivity to the weak rates is high but (within standard theory), the uncertainty is 
governed by the neutron lifetime which is now known with reasonable precision. Indeed, the 
relative uncertainty on r„ is of the order of 5xl0 -3 [4] affects 4 He abundance by (Table [1]) 
—0.73 x 5 x 10 -3 , a factor of ten lower than the observational uncertainty (Eq. [5|). 

The influence of the 1 H(n,7) 2 H rate was unexpected. The 7 Li final abundance depends 
strongly on the rate of this reaction while other isotopes are little affected. This unexpected effect 
can be traced to the increased neutron abundance at 7 Be formation time for a low 1 H(n,7) 2 H 
rate making its destruction by neutron capture, 7 Be(n,p) 7 Li(p,a) 4 He, more efficient (see Fig. 1 
in [22]). However, the few experimental informations available for this cross section at BBN 
energies are in good agreement with the calculations estimated to be reliable to within 1% 
error [20], 

The next most important reaction (Table [TJ is 3 He(a, 7) 7 Be as it is the path for the formation 
of 7 Li at high density. Hence, the 7 Li abundance is directly proportional to this rate, which has 
long been a subject of debate. Systematic differences in the measured cross section were found 
according to the experimental technique: prompt or activation measurements. Thanks to the 
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Figure 1. Astrophysical S-factor for the 3 He(a,7) 7 Be reaction, adapted from Fig. 3 in Di Leva 
et al. [24J. The data are from Refs. |23l [23] , the dot-dashed and solid curves are respectively the 
previously used Descouvemont et al. [21] (anterior to the displayed data [231 E|) or the adopted 
Cyburt & Davids [25] (including recent data [23]) fits. The dashed vertical lines represent the 
Gamow window at 1 GK. 



recent experimental efforts [23U24j . in particular at LUNA at the Laboratori Nazionali del Gran 
Sasso, the two methods provide now results in agreement, within each others error bars. With 
this new experimental data, Cyburt &: Davids [25] calculated the S-factor which is significantly 
higher than the Descouvemont et al. [21] R-matrix fit, done before these new data were available 
(Fig. [1]). This explains the higher 7 Li primordial abundance obtained in recent calculations. At 
high energy, the recent experimental data, in particular of Di Leva et al. [23], obtained by a 
third technique, the recoil mass separation, deviate from both fits. Theoretical explanations are 
available [26J, but this should not affect the S-factor at BBN energies. Nevertheless, one may 
note the relative scarcity of experimental data within the Gamow window (Fig. [TJ. 

The 2 H(d,n) 3 He reaction, also influential on 7 Li, was re-measured [together with 2 H(d,p) 3 H] 
by Leonard et al [27] . after the R-matrix analysis [21] was performed. The very precisely 
measured cross section is in perfect agreement with the R-matrix fit (Fig. [2]). 
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Figure 2. Experimental data for the, 2 H(d,n) 3 He reaction, S'-factor compared with the R- 
matrix fit of Descouvemont et al. |21j . (See this reference for reference to experimental data 
except for the more recent data from Leonard et al. [27J.) The dashed vertical lines represent 
the Gamow window at 1 GK 



2.3. BBN primordial abundances compared to observations 

Figure [3] shows the abundances of 4 He (mass fraction), D, 3 He and 7 Li (in number of atoms 
relative to H) as a function of the baryonic density. The thickness of the curves reflect the 
nuclear uncertainties. They were obtained by a Monte-Carlo calculation using for the nuclear 
rate uncertainties those obtained by [21] with the notable exception of 3 He(a,7) 7 Be |25j and 
1 H(n,7) 2 H |20| . The horizontal lines represent the limits on the 4 He, D and 7 Li primordial 
abundances deduced from spectroscopic observations. The vertical stripe represents the baryonic 
density deduced from CMB observations [5]. The concordance between BBN and observations 
is in perfect agreement for deuterium. Considering the large uncertainty associated with 4 He 
observations, the agreement with CMB+BBN is fair. The calculated 3 He value is close to its 
galactic value showing that its abundance has little changed during galactic chemical evolution. 
On the contrary, the 7 Li, CMB+BBN calculated abundance is significantly higher than the 
spectroscopic observations : from a factor of ~3[29j when using the Descouvemont et al. library 
|21| only and the Ryan et al. observations [30J (dotted lines in lower panel of Fig. |3}, to a 
factor of «5 [281 EI] when using the new rates and Li observations [7]. Table [2] displays the 
comparison between BBN abundances deduced from the WMAP results and the spectroscopic 



observations. The origin of this lithium discrepancy between CMB+BBN and spectroscopic 
observations remains an open question. Note also that with the new determination of the 
4 He primordial abundances [144 H5j . the agreement for 4 He becomes marginal. As shown on 
the figure, an increase of the rate of expansion of the universe during BBN (simulated by an 
additional effective neutrino family) would improve the situation. 



Table 2. Yields at WMAP baryonic density compared to observations. 



Cyburt el al 2008[28j CV10 [S] Observations Factor 



4 He 


0.2486±0.0002 


0.2476±0.0004 


0.2534±0.0083[I5] 


xl0° 


D/H 


2.49±0.17 


2.68±0.15 


3.02 ±0.23 H2] 


xl0~ 5 


3 He/H 


1.00±0.07 


1.05±0.04 


i.i±o.2H5] 


xl0~ 5 


7 Li/H 


c 24+ - 71 

°- Z4 -0.67 


5.14±0.50 


1.58±0.31[7] 


xlO- 10 



3. BBN production of heavier elements 

Even though the direct detection of primordial CNO isotopes seems highly unlikely with the 
present observational techniques at high redshift, it is important to better estimate their 
Standard Big Bang Nucleosynthesis production. Hydrogen burning in the first generation of 
stars (Pop III stars) proceeds through the slow p-p chains until enough carbon is produced 
(through the triple-alpha reaction) to activate the CNO cycle. The minimum value of the initial 
CNO abundance that would affect Pop III stellar evolution is estimated to be as low as 10 -13 
(in number of atoms relative to hydrogen, CNO/H) for the less massive ones [32J. This is only 
two orders of magnitude above the Standard Big Bang Nucleosynthesis CNO yield, using the 
current nuclear reaction rate evaluations of Iocco et al. [33 1. In addition, it has been shown 
that Pop III stars evolution is sensitive to the triple- alpha ( 12 C producing) reaction and can 
be used to constrain the possible variation of the fundamental constants [M]. This reaction 
rate is very sensitive to the position of the Hoyle state, which in turn is sensitive to the values 
of the fundamental constants. The same mechanism could also increase the amount of CNO 
( 12 C) produced in BBN. In the same context of the variations of the fundamental constants, 8 Be 
(which decays to two alpha particles within ~ 10~ 16 s) could become stable if these constants 
were only slightly different. At BBN time, this would possibly allow to bridge the "A=8 gap" 
and produce excess CNO [35J. To determine how significant would be this excess, one needs to 
know the standard BBN production of the CNO elements. 

The production of CNO isotopes has been studied in the context of standard and 
inhomogeneous BBN. The most relevant analysis comes from Iocco et al. [33J who included 
more than 100 nuclear reactions and predicted a CNO/H abundance ratio of approximately 
6 x 10~ 16 , with an upper limit of 10 -10 . 

The main difficulty in BBN calculations up to CNO is the extensive network needed, including 
n-, p-, a-, but also d-, t- and 3 He-induced reactions. Most of the corresponding cross sections 
cannot be extracted from experimental data only. This is especially true for radioactive tritium- 
induced reactions, or for those involving radioactive targets like e.g. 10 Be. For some reactions, 
experimental data, including spectroscopic data of the compound nuclei, are just nonexistent. 
Hence, for many reactions, one has to rely on theory to estimate the reaction rates. Previous 
studies lack documentation on the origin of the reaction rates, but have apparently extensively 
used old and unreliable prescriptions to estimate many of them. A detailed analysis of all 
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Figure 3. Abundances of 4 He (mass fraction), D, 3 He and 7 Li (by number relative to H) 
as a function of the baryon over photon ratio rj (or {l^-h 2 .) showing the effect of nuclear 
uncertainties [31J. The vertical stripe corresponds to the WMAP baryonic density [5] while 
the horizontal area represent the adopted primordial abundances (dotted lines those adopted in 
CV10 [31]). The dashed curves represent previous calculations [29] before the re-evaluation [25] 
of the 3 He(a,7) Be rate. The dot-dashed lines corresponds, to 4 effective neutrino families. 



reaction rates and associated uncertainties would be desirable but is impractical for a network 
of ~400 reactions. So we first performed a sensitivity study to identify the most important 
reactions, followed by dedicated re-evaluations. For this purpose, we used, as a first guess, the 
more reliable reaction rate estimates provided by the TALYS code |36j . 

3.1. Sensitivity study 

To extend our BBN network, we need the neutron, proton, deuterium, tritium, 3 He and a-particle 
capture cross sections on targets in the A=l to 20 range. As a first approximation, we used the 
results from TALYS for rates that are not available in the literature (the full list of references can 
be found in Ref. |37J). By comparing TALYS results with experimentally determined reaction 
rates [38l [39] , we observed that, even for very light elements like Li, TALYS globally provides 
predictions, "accurate" within 3 orders of magnitude, in the temperature range of interest here. 
Hence variations of these theoretical rates by up to three orders of magnitude can in a first 
step be used in our sensitivity analysis. Sensitivity studies have already been performed for the 
reactions involved in 4 He, D, 3 He and 7 Li production [40 [ I41 [ |29] but here, we will concentrate on 
the C, N and O isotopes production. To estimate the impact of the reaction rate uncertainties 
on Standard Big Bang Nucleosynthesis, we perform for each reaction six additional calculations, 
changing its rate by factors of 0.001, 0.01, 0.1, 10, 100 and 1000, and calculate the relative 
change in CNO abundances. (Mass fractions of isotopes with A>12 are added together into 
CNO.) Tables to O display, reactions for which the relative changes in 6 Li, 9 Be, n B and CNO 
abundances are larger than 20%. The last column contains the reference for the origin of the 
reaction rate used for the sensitivity study (the TALYS rates are available electronically [42]); 
i.e. before reaction rate re-evaluation (next section). 



Factors 


0.001 


0.01 


0.1 


10. 


100. 


1000. 


X < (TV > 


Reaction 




Fractional change in 


CNO 




Ref. 


Y Li(d,7) y Be 


1.00 


1.00 


1.00 


1.01 


1.11 


2.10 


TALYS 


7 Li(d,n)2a 


1.66 


1.65 


1.55 


0.28 


0.06 


0.02 


03] 


7 Li(t,n) 9 Be 


0.99 


0.99 


0.99 


1.10 


2.14 


11.7 


mi 


7 Li(t,2n)2a 


1.00 


1.00 


1.00 


0.99 


0.91 


0.53 


US] 


8 Li(n, 7 ) 9 Li 


1.00 


1.00 


1.00 


1.01 


1.06 


1.62 


02] 


8 Li(t,n) 10 Be 


1.00 


1.00 


1.00 


1.00 


1.02 


1.23 


TALYS 


8 Li(a, 7 ) 12 B 


1.00 


1.00 


1.00 


1.01 


1.11 


2.15 


TALYS 


8 Li(a,n) n B 


0.89 


0.89 


0.90 


1.97 


11.2 


78.1 


03 


9 Li(a,n) 12 B 


1.00 


1.00 


1.00 


1.01 


1.08 


1.73 


TALYS 


10 Be(a,n) 13 C 


1.00 


1.00 


1.00 


1.00 


1.03 


1.28 


TALYS 


n B(n, 7 ) 12 B 


0.91 


0.91 


0.92 


1.81 


9.91 


87.7 


06] 


n B(d,n) 12 C 


0.70 


0.71 


0.73 


3.67 


30.2 


280. 


TALYS 


n B(d,p) 12 B 


0.99 


0.99 


0.99 


1.08 


1.83 


9.33 


TALYS 


n B(t,n) 13 C 


1.00 


1.00 


1.00 


1.01 


1.12 


2.17 


TALYS 


n C(n, 7 ) 12 C 


1.00 


1.00 


1.00 


1.01 


1.08 


1.75 


06] 


n C(d,p) 12 C 


0.99 


0.99 


0.99 


1.05 


1.55 


5.67 


TALYS 


12 C(t,a) n B 


1.00 


1.00 


1.00 


1.00 


0.97 


0.75 


TALYS 


13 C(d,a) n B 


1.00 


1.00 


1.00 


0.96 


0.84 


0.75 


TALYS 



Table 3. Sensitivity of the most important reactions for CNO production in BBN, to reaction 
rate variations around initial test rates from references (last column). 



Factors 


0.001 


0.01 0.1 10. 100. 


1000. 


X < (TV > 


Reaction 




Fractional change in CNO 




Rcf. 


3 He(t,7) e Li 
4 He(d, 7 ) 6 Li 


1.00 
0.004 


1.00 1.00 1.03 1.31 
0.013 0.010 9.97 99.7 


4.11 
995. 


[48] 
[49] 



Table 4. Same as Table [3] but for 6 Li. 



Factors 


0.001 


0.01 


0.1 


10. 


100. 


1000. 


X < (TV > 


Reaction 




Fractional change in 


CNO 




Ref. 


Y Li(d,7) y Be 


0.83 


0.83 


0.85 


2.52 


17.7 


170 


TALYS 


7 Li(t,n) 9 Be 


0.52 


0.53 


0.57 


5.29 


48.2 


477. 


m 


7 Li( 3 He,p) 9 Be 


1.00 


1.00 


1.00 


1.04 


1.45 


5.49 


TALYS 


7 Be(d,p)2a 


1.01 


1.01 


1.01 


0.95 


0.67 


0.38 


m 


7 Be(t,p) 9 Be 


0.65 


0.65 


0.69 


4.15 


35.6 


345. 


TALYS 



Table 5. Same as Table [3] but for 9 Be. 



Factors 0.001 0.01 0.1 10. 100. 1000. x < av > 



Reaction 




Fractional change 


in CNO 


Ref. 


3 He(t,np) 4 He 


1.00 


1.00 1.00 1.00 


0.97 0.79 


[50] 


7 Be(d,p)2a 


1.01 


1.01 1.01 0.93 


0.55 0.11 


[50] 


n C(n,a)2a 


1.16 


1.16 1.15 0.40 


0.01 0.0001 


[46] 



Table 6. Same as Table but for n B. 



The examination of Table [3] shows that, among the ~400, only a few reactions have a strong 
impact on the CNO final abundance. The CNO production is significantly sensitive (more than 
by a factor of about 2) to several reaction rates. In particular, these include: 7 Li(d,n)2 4 He, 
7 Li(t,n) 9 Be, 8 Li(a,n) n B, n B(n,7) 12 C, n B(d,n) 12 C, n B(d,p) 12 B and n C(d,p) 12 C. The impact 
of Li(d,n)2 4 He is unexpected and should be compared to the influence of 1 H(n,7) 2 H on 7 Li 
(sec 22] ). Indeed, when increasing the 7 Li(d,n)2 He reaction rate by a factor of 1000, even 
though the 4 He, D, 3 He and Li final abundances are left unchanged, the peak 7 Li abundance at 
t Ri200 s is reduced by a factor of about 100 (see Fig. 15 in [37]), an evolution followed by 8 Li and 
CNO isotopes. From this table, we can deduce that the main nuclear paths to CNO (see also 
[33]) proceeds from the 7 Li(a,7) n B reaction followed by n B(p, 7 ) 12 C, n B(d,n) 12 C, n B(d,p) 12 B 
and n B(n,7) 12 B reactions. Another nucleosynthesis path starts with Li(n,7) 8 Li(a,n) 11 B. (Note 
that primordial n B is produced by a different path: the late decay of n C.) 

The examination of Tables H] to show that the most important reactions for 6 Li, 9 Be 
and n B productions are 4 He(d,7) 6 Li, 7 Li(t,n) 9 Be, 7 Be(t,p) 9 Be, 7 Li(d,7) 9 Be and n C(n,a)2a. 
However, some of them have been measured ( 4 He(d,7) 6 Li [49] and 7 Li(t,n) 9 Be |44[ 151] ) so 
that the uncertainties on their cross sections are small compared to the range explored in 
this sensitivity study. Considering the tiny production of these isotopes in Standard BBN 
(next section), compared to present day observational techniques, these uncertainties are not 
important. 

Our sensitivity study, with the extended network, also included 4 He, D, 3 He and 7 Li but 
no new important reaction was found besides those previously identified [291 HI], essentially the 
7 Be(d,p)2a reaction. Indeed, if this reaction rate were higher by a factor of ~100, 7 Li abundance 
would be brought down to the observed level |29| . An experiment, performed at Louvain-la- 
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Figure 4. Abundances in number of atoms relative to H for Li, Be, B, B and CNO 
isotopes. For CNO, the hatched area represent the present estimated uncertainty. For 6 Li, the 
dotted lines shows the uncertainty related to the 4 He(d,7) 6 Li reaction before the GSI experiment 
|49| (solid line )and the hatched area shows the most significant 6 Li detection in HD 84937 star 
[6Q1E]. 



Neuve did not find such an enhancement in the cross-section integrated over the Gamow window 
|52j . Afterwards, Cyburt & Pospelov [53] proposed a resonance enhancement of the cross section 
that could have been left undetected by this experiment. Later, a dedicated experiment at Oak 
Ridge [53] did not find such a resonance, in the 7 Be+d channel. Then, Kirsebom & Davids [55 1 
pointed out that the properties of the corresponding 9 B level had been measured [56]. When used 
in the reaction rate and subsequent BBN calculation, the 7 Li depletion was found insignificant 
|55j (<4%). The 7 Be+ 3 He channel was found promising by Chakraborty et al. [57J since the 
spectroscopy of the compound nucleus 10 C is deficient in the Gamow window. But as in the 
7 Be+d channel, the required level properties are at the fringe of standard nuclear physics, as 
shown by Broggini et al. [58J. In addition, "missing" 10 C levels, were not found in a dedicated 
experiment, recently carried out in Orsay [59] . It seems, then that the possibility of a nuclear 
solution to the lithium problem, has been ruled out. 

3.2. Results 

For the few reactions that were identified to have an impact on Standard Big Bang 
Nucleosynthesis, we (re-)evaluated the rates issued from TALYS, or other sources whose 
references are listed in the last columns of Table [3] to Table [6j We were able to collect sufficient 
experimental data, or theoretical constraints, to derive new reaction rates with associated 
uncertainties, much reduced with respect to our initial three orders of magnitude variation. 
In some cases these new rates differ from the previous ones by large factors but changes 
compensate each other (e.g. 11 B(d,n) 12 C and 11 B(d,p) 12 B). In the meantime, the important 
8 Li(a,n) 11 B reaction rate, was independently re-evaluated by La Cognata & Del Zoppo |61j . but 
it affects CNO production by less than 2/^, when compared to our own re-evaluation. We hence 
confirm that the CNO Standard Big Bang Nucleosynthesis production is CNO/H ~ 0.7 x 10 -15 
(number of atoms relative to H). Our present analysis does not allow us to precisely quantify 
the uncertainty on this result, but from the inspection of Tabled and assuming a factor of ten 
uncertainty on the re-evaluated reaction rates, we can estimate the range of CNO/H values to 
(0.5 — 3.) x 10~ 15 . These results are consistent with those of Iocco et al. |33j but slightly lower 
than those of Vonlanthen et al. [62]. Detailed CNO, B and Be isotopic abundances, compared 
with those from Iocco et al. [33] can be found in Tabled together with H, He and Li isotopic 
abundances compared with our previous work. As expected, the extension of the network does 
not alleviate the Li discrepancy between calculations and observations. No uncertainty is 
available yet for the extended network result but, for 4 He, D, 3 He and Li, they are within the 
reduced network uncertainties |31| . The differences in the central values are essentially caused 
by the small evolution of the baryonic density, following the progress in WMAP data reduction 
over the years. 

For 6 Li, now that an increased 2 H(a, 7) 6 Li cross section is excluded[49], the BBN 6 Li yield 
( 6 Li/H~ 10 -14 ) at WMAP baryonic density is about two orders of magnitude below the reported 
observations in some halo stars, that nevertheless have to be confirmed. 

The CNO Standard Big Bang Nucleosynthesis production is found to be little sensitive 
to the baryonic density of the Universe as shown on Figure S where the 6 Li, 9 Be, 10 B and 
n B abundances, results of our calculations, are also depicted. Even when considering our 
estimated uncertainty, the primordial CNO abundance is too low to have an impact on Pop III 
stellar evolution. Those first stars having high masses and consequently short lifetimes are not 
observed today. The high C and O abundances observed in extremely low metallicity Pop II 
stars are expected to have been synthesized by massive Pop III stars making the determination 
of primordial CNO abundances from observations presently out of reach. 



2 1.4% indeed, and not a factor of 1.4 (a typo in [37] ~). 



Table 7. Primordial abundances at (slightly evolving) WMAP baryonic density, with an 
extended network and after reaction rate re-evaluations, compared to previous works. 





CV10 [31 


CGXSV [37J 




0.2476±0.0004 


0.2476 


D/H (xl(T 5 ) 


2.68 ± 0.15 


2.59 


3 He/H (xl(T 5 ) 


1.05±0.04 


1.04 


7 Li/H (xlCT 10 ) 


5.14±0.50 


5.24 


6 Li/H (xlO~ 14 ) 


1.309] 


1.23 




Iocco et al. [33 


CGXSV [37] 


y Be/H (xlO" iy ) 


2.5 


9.60 


10 B/H (xnr 21 ) 




3.00 


11 B/H (xl0~ 16 ) 


3.9 


3.05 


12 C/H (xlO -16 ) 


4.6 


5.34 


13 C/H (xl0~ 16 ) 


0.90 


1.41 


14 C/H (xHT 21 ) 


13000. 


1.62 


14 N/H (xlO- 17 ) 


3.7 


6.76 


15 N/H (xlO" 20 ) 




2.25 


16 0/H (xlO^ 20 ) 


2.7 


9.13 


CNO/H (xlO" 16 ) 


6.00 


7.43 



4. Conclusions 

The baryonic density of the Universe as determined by the analysis of the CMB anisotropies is 
in very good agreement with Standard BBN compared to D primordial abundance deduced from 
cosmological cloud observations. However, it disagrees with lithium observations in halo stars 
by a factor that has increased with the availability of improved nuclear data and astronomical 
observations. Presently, the favored explanation is lithium stellar depletion, but the larger 
needed depletion factor is hardly compatible whit the thin observed plateau. It is hence essential 
to determine precisely the absolute cross sections important for 7 Li nucleosymthesis (Table H}. 

Nevertheless, primordial nucleosynthesis remains an invaluable tool for probing the physics 
of the early Universe. When we look back in time, it is the ultimate process for which we, 
a priori, know all the physics involved. Hence, departure from its predictions provide hints 
for new physics or astrophysics. Gravity could differ from its general relativistic description, 
for instance a scalar field, in addition to the tensor field of general relativity (GR), appears 
naturally in superstring theories. That would affect the rate of expansion of the universe and 
hence BBN (see [65j and Ref. [63] for a review) . Coupled variation of the fundamental couplings 
is also motivated by superstring theories (see Ref. [65] for a review) . However, the impact of these 
variations on the nuclear reaction rates is difficult to estimate, as in general, nuclear physics uses 
phenomenological models, whose parameters are not explicitly linked to fundamental constants. 
The decay of a massive particle during or after BBN could affect the light element abundances 
and potentially lower the 7 Li abundance (see e.g. [66]). Negatively charged relic particle, like 
the super symmetric partner of the tau lepton, could form bound states with nuclei, lowering 
the Coulomb barrier and hence leading to the catalysis of nuclear reactions (see e.g. [67\ 68J). 
Annihilation of dark matter during BBN, e.g. injecting extra neutrons, could also modify the 
primordial abundances [69l [70] . 

We have extended our network up to the CNO region and performed a sensitivity study to 
identify the few reactions that could affect the A>7 isotope yields and re-evaluated their rates. 
The CNO isotope production was found to be in the range CNO/H = (0.5 — 3.) x 10~ 15 , not 
sufficient to have an impact on the evolution of the first stars. It is nevertheless a reference 



value for comparison with non-Standard Big Bang Nucleosynthesis CNO production e.g. in the 
context of varying constants. In this particular case, even with a faster triple-alpha reaction 
rate or a stable 8 Be, the C(NO) production remains ~6 order of magnitude [37] lower than the 
Standard Big Bang Nucleosynthesis value reported here. 

Last but not least, we stress here the importance of sensitivity studies in nuclear astrophysics 
that have been done, e.g. in the context of novae [71], X-ray burst [72] or massive stars [73]. Even 
in the simpler context of BBN without the complexity (e.g. mixing) of stellar nucleosynthesis, 
it would have been very unlikely to predict the influence of the 1 H(n,7) 2 H reaction on 7 Li nor 
of the 7 Li(d,n)2 4 He reaction on CNO. 
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